Level correlations driven by weak localization in 2-d systems 
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We consider the two-level correlation function in two-dimensional disordered systems. In the non- 
ergodic diffusive regime, at energy e > Ec {Ec is the Thouless energy), it is shown to be completely 
determined by the weak localization effects, thus being extremely sensitive to time-reversal and spin 
symmetry breaking: it decreases drastically in the presence of magnetic field or magnetic impurities 
and changes its sign in the presence of a spin-orbit interaction. In contrast to this, the variance of 
the levels number fiuctuations is shown to be almost unaffected by the weak localization effects. 
(Submitted to Phys. Rev. Lett. 19 Septmeber 1994) 
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A remarkable observation of Wigner and Dyson [1] was 
the existence of universal statistics that govern energy 
level spectra in a wide variety of quantum systems [2]. 
The Wigner-Dyson statistics of eigenvalues in the invari- 
ant ensembles of random matrices are applicable also to 
disordered electron systems in a metallic phase [3-6]. The 
level repulsion, characteristic of these statistics, results 
from the overlaping of one-electron states corresponding 
to different energies. With the disorder increasing such 
systems undergo the Anderson metal-insulator transition 
[7]. On an insulator side of the transition, uncorrelated 
energy levels are described by the Poisson statistics. 

Recently, the existence of the third universal spec- 
tral statistics, applicable to the disordered systems in a 
critical regime near the Anderson transition point, has 
been predicted [8-11]. These statistics emerge [10] in the 
crossover region near the transition, where electron diffu- 
sion is anomalous, and turn out to be completely new [9] 
and different from the Poisson and Wigner-Dyson ones, 
rather than a universal hybrid of both, as has been earlier 
conjectured [8]. They are still characterized by the level 
repulsion [11], albeit it is weaker than in the metallic 
region. It is essential that nontrivial level correlations 
[9,10] turn out to be driven by small scale-dependent 
corrections to the conductance which by itself is almost 
scale-independent near the Anderson transition. 

The conductance has only a weak scale-dependence 
also in the case of a two-dimensional disordered system 
in a weak-localization (WL) regime. There is no metal- 
insulator transition at d = 2, and, with the size increas- 
ing, such a system crosses over to a strong-localization 
regime [12]. So one can expect some analogies between 
the level correlations in this case and those in the criti- 
cal regime in the higher dimensionality. However, while 
the one-parameter scaling hypothesis [12] was taken for 
granted to make considerations in the critical regime pos- 
sible [9,10], the WL regime at d = 2 may be treated within 
a rigorous and conjecture-free perturbative approach. 

In this Letter we consider level correlations in a 2d 
disordered electron system and show that they are to- 
tally governed by the WL corrections to the conductance, 
when the distance between the levels, e, is in the region 



Ec<e< h/gr (1) 

This is another of a few known examples (as an anoma- 
lous magnetoresistance or the Aharonov-Bohm effect [7]) 
where the weak-localization corrections determine the 
main effect. In Eq. (1), Ec = HDq/Lq is the Thouless 
energy, Lq is the sample size, Dq = tVpl2 is the diffusion 
coefficient, 1 = ii^r is the elastic scattering length, and 
g is the conductance in units of e'^/Tr/i. 

We consider the two-level correlation function: 
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Here (. . .) denote averaging over all the realizations, and 

p{E)=pL-^Y.^^^-^^) (3) 

n 

is an exact density of states for a particular realization of 
disorder, where e„ are exact energy levels, and p is a spin 
degeneracy. For E, + £ in the bulk of conduction band, 
{p[E)) = p is energy-independent, and Re{s) = is 
a function of (e) only. 

The averaging in Eq.(2) is performed within the impu- 
rity diagrammatic technique [13] with p{E) represented 
by [ip/ttL'^) JlmG'^[r,r; E)dr (here are exact re- 
tarded or advanced one-particle Green's functions). It 
is convenient to use a representation in which slow dif- 
fusion modes are explicitly separated from fast "ballis- 
tic" ones [14,15]. In Fig. 1, there are shown the lowest 
order (one-loop) diagrams containing one or two diffu- 
sion or cooperon propagators (shown by wavy lines which 
correspond to the ladder series in the conventional tech- 
nique [5]). Triangular and square " Hikami boxes" repre- 
sent the motion at the ballistic scale as all the averaged 
Green's functions G [r—r',£) (represented by edges of 
the boxes) decay as e"''^ l'^^. At the diffusion scale, 
these boxes are reduced to certain constants (see Ref. 
[16] for review). Taking into account the total number 
of cooperon and diffuson contributions that depends on 
the universality class, one obtains for the one- and two- 
diffuson diagrams: 
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FIG. 1. One-loop diagrams contributed to J?(e) 



1 



A2 



1 



(4a) 
(4b) 



oyj, and nx,y run over all the in- 



where q= [2^ / Lo){nx , riy 
tegers, /3= 1, 2, or 4 is for the Dyson orthogonal, unitary, 
and symplectic ensembles, respectively, A=p/pL'^ is the 
mean level spacing. 

The two-diffuson diagram of Fig. lb is well-known to 
dominate at c! > 2 in the interval A ^ e ^ Ti/r. Indeed, 
in the ergodic region, £ <^ Ec, where all the terms with 
q ^ may be neglected in Eq. (4b), it gives 
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which is the envelope of the Wigner-Dyson correlation 
function. The contribution of the one-diffuson diagram, 
Eq. (4a), is of order {1/ LoYilkp)'^^'^'^'' in this region, 
i.e. totally negligible. In the diffusive region, Ec <^ £ <^ 
h/r, the summation in Eqs. (4) may be changed by in- 
tegration, and the two-diffuson diagram, Eq. (4b), gives 
the well-known result [5]: 
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Here the numerical coefficient Cd depends only on d. 

For d > 2, Eq. (6) makes the main contribution into 
i?(e), as that of the diagram of Fig. la turns out to be 
smaller by (Z/i^)^"'^, where Ll = hDo/£. It becomes rel- 
evant for the correlation function (2) only at the ballistic 
energy scale e ~ h/r, where the appropriate diffusion 
length Ls is of order i and the diffusion approximation 
is no longer valid. 

For d = 2 such an estimation fails as C2 = in Eq.(6). 
A similar cancellation (of a contrubution that would be 
main if not a vanishing coefficient) happens at the mo- 
bility edge for d > 2, although much more subtle con- 
siderations were required [9]. On the face of it, what is 
left at c! = 2 is the contribution of the diagram of Fig. 
la and the next order correction for the diagram of Fig. 



lb (arisen when the corrections to the boxes oc {qi)^ are 
taken into account). It gives: 
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This result corresponds to that obtained by Altland and 
Gefen [17] who have thoroughly considered the level cor- 
relations in the ballistic regime (in their technique both 
the diagrams of Fig. 1 were described by a single expres- 
sion so that no explicit cancellation has been seen). In 
the diffusive regime, it gives only a logarithmic depen- 
dence on £ with a small amplitude, rA/gh ~ (l/kpLoY . 

A surprising situation at c! = 2 is that in a wide re- 
gion, Eq. (1), the contribution of the one-loop diagrams, 
Eq. (7), is small compared to that of the two-loop dia- 
gram containing four propagators (Fig. 2a). The shaded 
Hikami box is oc [15], so that this diagram describes 
the WL correction to the diffusion coefficient Dq, in Eqs. 
(4). It contains a cooperon, thus being absent in the (3 = 2 
case where the next-order (three-loop) diagrams should 
be considered. Therefore, to find the WL contribution to 
the correlation function (2), we calculate the correction 
to Do, given for /3 = 1, 4 by 
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and substitute D = Dq +SD into the integral in Eq. (6). 
Naturally, its real part would not change the zero result 
at d = 2. So it is the imaginary part of the WL correction 
that governs the level correlations in the region (1). With 
the same accuracy as in Eq. (6), 
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FIG. 2. Two-loop diagrams contributed to J?(e). The dia- 
grams (b) and (c), that describe renormalization of the effec- 
tive vertex, cancel each other. The diagram (a) yields the WL 
correction to the diffusion coefficient. 
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Here 77(5, /3) = (/3-2)-i for /3 = 1, 4, and = -l/2g for 
/3 = 2 where the above result is obtained with due regard 
to the third-loop corrections, and thus contained an ad- 
ditional factor of 1/g. Calculating the integral (6) with 
allowance for the correction f9), we obtain 
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There exist two more two-loop diagrams, Fig. 2b, c, but 
they exactly cancel each other. It has been expected, 
as these diagrams would yield a singular correction to a 
renormalized triangular vertex of the diagram of Fig. lb. 
The vertex has zero scaling dimensionality, and should 
remain unrenormalized, as its renormalization would vi- 
olate the particle conservation law. Thus all the vertex 
corrections must cancel each other, as has been noted in 
Ref. [9] for the diagrams at the mobility edge. The can- 
cellation of the diagrams in Fig. 2b, c provides a pertur- 
bative proof of the absence of the vertex renormalization. 

In the region (1) this contribution dominates over that 
of Eq. (7), provided that g < h/EcT ~ {Lo/lf. (For 
/3 = 2, it dominates for ^ h/g^T, such region existing 
for g < Lq/I). This condition may be violated for a very 
small or a very pure sample which is almost ballistic. In 
general, however, it is easily satisfied. The two contribu- 
tions are matching for e ~ hrj/gr, and the latter becomes 
dominating for the higher energies, in the 'quasi-ballistic' 
region hrj/gr ^ e ^ h/r. There is, however, a paramet- 
ric mismatching of the results of Eqs. (5) and (10) at 
the boundary between the ergodic and diffusive regime, 
e ~ Ec, where the ergodic expression, Eq. (5), is of order 
— 1/g^ while diffusive one is of order —1/g^ for (3=1, 
-1/g^ for /3 = 2, and +l/g^ for /3 = 4. The reason is 
that one can justify neglecting all the g terms in 
the sum (4b) only for £ <^ Ec, and changing summa- 
tion by integration only for e ^ Ec- There exists some 
transient region near e ^ Ec where -R(e) crosses over 
from the Wigner-Dyson behavior, Eq. (5), to the WL 
behavior, Eq. (10). In this region one should evaluate 
accurately the sum in Eq. (4b). Changing there to the 
dimensionless variables n and s, and using the transfor- 
mation = Jq ae~°'^da, one performs the summation 
over rix and Uy in Eq. (4b) to obtain the result in terms 
of the elliptic theta-function [18] ^3: 
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The asymptotical values of (p are given by 




2 exp(-7r^/a) 



a < 1; 
a > 1. 
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For s <^ g, Eq. (11a) is thus reduced to the ergodic result 
(5), while for s ^ g, Re/(s) oc exp[— 7r-^2s/g]. There- 
fore, for large s using the accurate summation instead 
of integration gives only the exponentially small correc- 
tion to the WL result of Eq. (10). (Note that with the 
accuracy up to the same exponent, Im/(s) = tt/s, so 
that using the integration for calculating the WL correc- 
tion is justified). In the transient region of the width 

= 2Ecln^ g, the contribution of Re/(s) is important 
(and dominant for |e — <C £„), and it describes a 
smooth and universal crossover from the Wigner-Dyson 
to the WL regime in R[s). 

Now we can combine the results in the ergodic re- 
gion, Eq. (5), the WL region, and the transient region 
described above to represent them as 
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Only for s > [h/grA), this contribution becomes smaller 
than the quasi-ballistic one, Eq. (7), found in Ref. [17]. 

The two-level correlation function R{s) studied above 
determines the variance E2 = {{N — N)'^) of the number 
of levels N in an energy strip of the given width E = NA 
(here (N) = N). Its derivative is given by 
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(14) 



In order to calculate the r.h.s of Eq.(14), it is convenient 
to use the representation of R{s) in the form of a sum, as 
in Eq. (lib), and integrate first over s. Then, neglecting 
the WL corrections, one obtains for N ^ g: 
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so that E2 ~ N/g [5]. Allowing for the WL term, one 
finds with the logarithmic accuracy: 
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As in the WL regime g~^hi[Loll) << 1, the variance 
E2 is almost linear in N in the whole region Ec < E = 
NA < h/r. Indeed, the WL logarithmic term in Eq.(16) 
could be important only for N > ge^ . For E ^ h/r, this 
gives ^ h/EcT ~ [Lq/I)^ while the WL approach is 
applicable in the opposite limit, g~^\n[Lo/l) ^ 1. 

Note that the AT-independent result of Eq. (15) is only 
valid in the diffusive region h/[TA) N ^ g. For 
N ^ 00 the integral in the r.h.s. of Eq. (14) must vanish. 

The reason is that there exists a sum rule which follows 
from the conservation of the total number of states Af. 
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Multiplying Eq. (2) by {p{E + e)) and integrating over 
all e with the identity L'^ J p[E + e) de = J\f , one obtains: 

//. + 00 
RE{e){p{E + e)) de = / R{s)ds = Q, (17) 
J — oo 

where the second identity in Eq.(17) results from (p) be- 
ing a constant in the relevant interval, e <C • 

By the derivation, the sum rule Eq.(17) holds exactly 
as long as the total number of states M is finite. How- 
ever, it may be violated if the thermodynamic (TD) limit, 
M ^ oo, were taken before the integration in Eq.(17). 

To illustrate this point let us consider an ensemble of 
diagonal random M x A/'-matrices. Suppose that ma- 
trix elements, independent random variables each 
having the same probability distribution: P{x) = 
for \x\ < M /2 and P{x) = otherwise. Using the defini- 
tions (2, 3), one immediately finds the correlation func- 
tion, i?„(s) = S{s)-Af-^e{Af/2-\s\), where e{x) is a step 
function. For finite Af, -Ro(s) obeys the sum rule (17), 
since a positive contribution of the J-function is cancelled 
by a negative one given by the flat small "tail" in -R(,(s). 
The universal Poisson distribution, Ro{s) = S[s), arises 
only in the TD limit and obviously does not obey the 
sum rule. 

The above toy model illustrates what happens in the 
present problem. The integral in Eq. (14) is really van- 
ishing in the limit N oo, as required by the sum rule, 
due to the negative contribution, [—f3g)~^, of the long 
and small "ballistic" tail in R{s) given by Eq. (7) for 
er/h ~ 1. If the TD limit, Lq oo, is taken in R{s) this 
ballistic tail vanishes, and the sum rule breaks down. 

A situation similar to the present 2d case takes place 
at c! > 2 at the mobility edge where the sum rule for the 
limiting function R{s) breaks down for the same reason. 
Thus, in contrast to the statement made in Ref. [9], it 
could not be applied to calculating the variance. This 
leads to a linear in N term in the variance [19], missed in 
Ref. [9], which is similar to that in Eq. (16). A difference 
is that g ~ 1 at the mobility edge so that the coefficient 
of proportionality is just a certain number. The long- 
range levels correlations lead to the contribution to the 
variance oc JV^ ( 7 is a universal critical exponent), so 
that E2 = AN -\- BN'^ where A and B are contributed by 
all the diagrams and thus cannot be exactly determined. 

We conclude that in the diffusive region of Eq. (1) the 
level correlations, Eq. (10), are totally governed by the 
weak localization effects while the level number variance, 
Eq. (16), is governed by the total conductance with the 
WL effects resulting only in small corrections. 
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